Wormlike chain or tense string? A question of resolution 
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It is shown that a wormUke chain, i.e., a filament with a fixed contour- length S and a bend- 
ing elasticity k, attached to a frame of length L, can be described — at low resolutions — by 
the same type of elastic free-energy as a tense string. The corresponding tension is calcu- 
lated as a function of temperature, L, k and S. 



Elastic free-energies are used to describe the distorted states and fluctuations of soft- 
matter ( Ctiaikin and Lubensky, 1995| ). The latter includes liquid-crystals ( |de Gennes and Prost 



1993 ), surfactant membranes or interfaces ( [Israelachvili, 1991 ), and polymers ( Doi and Edwards 



1989). In this paper, we investigate how the elastic free-energy functional of a filament with a fixed 



contour-length and a bending rigidity (wormlike chain) depends on the spatial resolution at which 
it is described. We perform a finite renormalization group (RG) iteration in order to determine the 
free-energy associated with the length of the filament that is hidden at low-resolution. A similar 
study, focused on surfactant membranes, was recently published in a short paper ( [Fournier et al 



20011) . 



The filaments considered here are known as wormlike chains ( [Kratky and Porod, 1949 ), or 
semi-flexible polymers ( poi and Edwards, 1989 ; Landau and Lifshitz, 1971). At temperature T, 
a filament with a bending rigidity k has a persistence length Lp = K/k^T ( [Landau and Lifshitz 



1971): for contour-lengths S ^ Lp it resembles a flexible rod, while for S" ^ L„ it resembles a 



wandering polymer. There are various natural realizations of semi-flexible polymers: actin cy 
toskeletal filaments, Lp ~ 4 — 17 /im ( [Kas et al., 199'| ), microtubules, Lp ~ 6mm ( [Gittes et al. 
1993[ ), and DNA, Lp ~ 50 - 350 nm pagerman, 1988[ ; [Marko and Siggia, 199"5[ ). Here we shall 



consider fluctuating filaments with S ^ Lp, however they will be stretched by attachment to a 
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frame of width L < 5. 

The paper is organized as follows: in Sec. I we recall how a coarse- graining procedure allows to 
define different free-energies associated with different (spatial) resolutions of the same system; in 
Sec. H we describe the high-resolution free-energy of a wormlike chain; in Sec. Ill we calculate the 
corresponding low-resolution free-energy, and we show that it coincides with that of a tense string; 



in Sec. IV we summarize and discuss our results. 



I. CONSTRUCTING ELASTIC FREE-ENERGIES 

The free-energy of a system in thermodynamic equilibrium with a heath bath at temperature T 
is defined by 

F[S] = -kBT In (y.^'^""')^ (1) 
\ies ) 

where /? = Xlk^T is the inverse temperature and i labels the microscopic states (microstates) of 
the system, Jii being the energy, or Hamiltonian of state i. The above sum Yl,%^s exp(— /^Wi), also 
called restricted partition function Z5, runs over different subsets S of the phase space F of all the 
microstates of the system. Each subset S defines a macrostate: it is a collection of microstates that 
we wish to consider macroscopically as a whole. F[S] is the free-energy of the macrostate S; it is 
either a functional or a function depending on the explicit form of S. For instance, in the standard 
thermodynamics of a fluid, the macrostates S are defined by specifying the number of particles N 
and the volume V of the fluid. Hence, F is a function of the variables T, V, N. 

It is most convenient to divide the phase space T into a collection of disjoint macrostates S. 
Then, since we know from statistical mechanics that the probability for the system to be in an 
arbitrary microstate i is given by 

P^ = ^e-^'^', (2) 

where Zr = Z^ier 6xp(— /JWj), the quantity 

P[S] = -le-'^^I'^l (3) 

is the probability to find the system in the macrostate S. Indeed, from Eq. ([l|) we have 
Z^^ exp(— /3F[5]) = J2ieS exp(— /J'Hj). In other words, F[S] plays the role of a Hamiltonian for 
the macrostates S. If the macrostates are well-chosen and large enough, the system will fluctuate 
little around the macrostate S* that minimizes F[S], i.e., around the most probable macrostate. 
The latter determines then the macroscopic "equilibrium state" of the system. 
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In general, for a soft-matter system, one chooses for the macrostates S the profiles of the order- 
parameter field: the latter is defined as the local average of some interesting material property, 
e.g., for a nematic^, Qij{r) = {mi{:x.)mj{x.) — ^Sij), where m(x) is the orientation of the molecule 
situated at point x, the average running over a small volume centered on r. Note that the size Aq ^ 
of this volume may be chosen at will and sets a cutoff Aq for the order-parameter field. In the case 
of a filament, one may simply choose the positions r(s) of its constituting elements. 

It is essential to note that we have a complete freedom to choose the "resolution" at which 
the phase space F is coarse-grained into the macrostates S. For example, on the one hand, we 
may be interested in the microscopic fluctuations of a bilayer membrane, as resolved at the scale 
of the membrane thickness (~ 40A). We would then define the free-energy F[S], as in Eq. (|l]), by 
summing over all the microstates i compatible with a given shape S of the bilayer. Specifying one 
of these microstates would require specifying all the coordinates, velocities and conformations of 
all the lipids within the bilayer. On the other hand, we may be micro-manipulating a membrane 
under an optical microscope. In this case, we would only be interested in the Fourier components 
of the membrane shape with wavevector q less than the typical wavevector of the visible light. We 
would then construct a free-energy functional F[S] by summing over all the microstates compatible 
with the poorly resolved shape S of the membrane ( [Fournier et ai, 2001 ). This would give us a 



tool more adapted to the experiment, since exp(— is directly the probability of observing S 
(up to a normalization constant). 

The free energy F[T,] = —k^T In J2ii:Y;^^Pi~P'^i) associated with the coarse-grained states S 
can be calculated also from 



F[S] = -A:Brin('^e-^^['5A 
\5eE / 



(4) 



since 



ies SGSie5 <sgs 

Equations (^) and (^) define a "coarse-graining" procedure that can be carried-out for any partition 
S of the phase space that is obtained by grouping the elements of a sub-partition S. Note that 
the RG consists in applying this procedure repeatedly while rescaling the lengths and the fields at 
each step ( poldenfeld, 199^ ). 



^ a fluid with a long-range non-polar orientational order of its constituent molecules. 
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It is usually impossible to directly perform the microscopic sum defining However, as- 

suming on the grounds of symmetry^ a simple form for F[S], it may be possible to determine 
F[T,] analytically from Eq. (^). The functional F[E,] will in general involve renormalized elastic 
constants ( phaikin and Lubensky, 1995 ) depending both on the elastic constants of F[S] and on 



temperature T (these new elastic constants take into account the effects of the thermal fluctuation 
occurring at shorter scales). However, may have a functional form fairly different from that 
of F[S], as we shall see in the following. 

It should be noted that this procedure, although akin to that of the RG, is different from 
the point of view of its scope. The scope of the RG is to identify the phases and critical points 
of a system by coarse-graining/rescaling infinitely. Here, the scope is to construct a free-energy 
functional adapted to the resolution at which the system is actually observed, or described (i.e., a 
tool fit to determine the shape and fluctuations of the system at this resolution). 

II. A WORMLIKE CHAIN AT HIGH RESOLUTION 

Let us consider a filament of fixed contour-length S, attached at its extremities to two fixed 
points, separated by a distance L along the x-axis (Fig. ||). We choose as macrostates S (see Sec. I) 
the microscopic shape of the filament, measured at the resolution that corresponds to its thickness. 
We restrict our attention to small deformations of the filament, hence we parameterize its shape by 
a normal displacement field r(x), with r -x = and x G [0, L]. Each macrostate S thus corresponds 
to a field r(2;). Let us introduce its Fourier transform Vq-. extending r(x) periodically in such a way 
that r(x + L) = r(x) for each x, we may write 

r(^) = 4FEr.e"'^ (6) 

Vi> g 

L 

dxr(x)e-^'?^, (7) 



1 

VL Jo 

with q = n X lir/L, where n is a positive or negative integer. Since the highest resolution for r(x) 
corresponds to the filament's thickness ~a, we require 

|g| < Ao = — . 8 

a 

The filament is assumed to be made of a fixed number of incompressible segments, hence its 



^ One usually represents the macrostates 5 by a small deformation field 0, and one writes F as the integral of a 
free-energy density expanded in power series of (p and its gradients. 



FIG. 1 Parameterization of the filament's shape. The disc indicates the plane normal to x, in which v{x) 
lives. In practice, we shall consider filaments that remain in the vicinity of the x-axis. 



contour-length 



is fixed. Here we have used the exact expression for the elementary arc-length ds = 
dxyjl + {dv/dxY and we have expanded the result up to second-order in y{x). In the following, 
assuming small tranverse excursions of the filament, we shall systematically expand all expressions 
in power-series of r(x), and retain only second-order terms. 

Using Eq. (^) and J^dx ex.p{iqx) = L6qfi, where 6 is the Kronecker symbol, we have 

S-L^ l^'l^'^l"- (10) 

|g|<Ao 

We have used the property r„g = r^* that follows from the fact that r(x) is real-valued. 

We assume that the free-energy of the filament, obtained in principle from the microscopic ma- 
terial properties through the procedure embodied in Eq. (||), is simply proportional to the integral 
of the square of the filament's curvature ( parris and Hearst, 1966 ; Kratky and Porod, 194g| ). More 
precisely, 

F[r] = oo otherwise. (12) 

Let us calculate the squared curvature up to second-order in r. From dM = dxx + dr, where M is a 
vector the extremity of which runs along the filament, we deduce the tangent vector t = dM/ds ~ 
[1 — ^{dr/dx)'^]ii. + dr/dx. Hence, the curvature is dt/ds ~ — [(d^r/dx^) • (dr/dx)]'x. + d'^r/dx^ and 
(dt/(is)2 ~ (d^r/dx^f. The functional in Eq. (Ill]) becomes, up to second order in r. 



F[r]^ tdx'-J^y = i^g^lr/. (13) 

^ / kl<Ao 
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III. THE SAME WORMLIKE CHAIN AT LOW RESOLUTION 

Let us now assume that we are not interested in the Fourier components of the filament's shape 
in the range A < \q\ < Aq (with A > 0). For instance, we may be observing the filament with an 
experimental setup of resolution 2ir/A. The coarse-grained macrostates S of interest (see Sec. I) 
correspond then in real space to 

R(x)= Y: ^,e'^\ (14) 

kl<A 

The field R(x), i.e., the poorly resolved shape of the wormlike chain, has a cutoff A lower than that 
of r{x). We are interested in the free-energy F\R\ defined as in Eq. (^) and not in F[r\. Indeed, 
exp{— /9F[R]} is, up to a multiplicative constant, the probability for the occurrence of a given 
low-resolution shape R(x) — with any possible wrinkle in the either inaccessible or uninteresting 
range A < |g| < Aq; on the contrary, exp{— /3F[R]} would be the probability for the occurrence of 
a given low-resolution shape R(x) having no high-wavevector wrinkles at all (which is useless since 
it corresponds to a very particular and unlikely event). Note that in a purely harmonic system — 
which is not the case here because of the fixed contour-length constraint — the probability weight 
of the wrinkles is actually independent of the low-resolution shape R(x); in this case -F'[R] = -F'iR]. 
Following the procedure embodied in Eq. (^), we have 

F[R] = -fcBTln /■ n dv'^dv';^5[L+ J2 ^^V.l' " ^ ) exp I -/? ^ l^q^lrA . (15) 

A«ir<Ao V |g|<Ao / V |9|<Ao / 

Here, and are the real and imaginary parts of r^, respectively. The integration variables [the 
wrinkles hidden in R(x)] span the range A < g < Aq of positive wavevectors only, because r_g and 
Tq are not independent (they are complex conjugates). In principle, the above calculation is valid 
only for small values of r^, however, thanks to the Gaussian character of the integrand, we take 
the integrals as running over ] — cxd, oo[. The Dirac delta- function takes into account the constraint 
on the filament's contour length [see Eq. ([T^)]. 

Let us call S the contour-length of the low-resolution shape R(x), i.e., the apparent contour- 
length. The latter, which is no longer fixed, is given by 

S^L+Y: ll'lr.l'. (16) 

kl<A 

Using this relation and the integral representation 6{x) = J^^dXexp{iXx) of the delta-function. 
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we may rewrite Eq. (15) in the form 

dX 



exp{-/3F[R]} 



n « 



A<(j<Ao 



exp 



A<|g|<Ao 



kl<Ao 



(17) 



Since = Rg for |g| < A, we have 



1 



exp{-/3F[R]} = exp | -/? ^ -ac/|R, 

kl<A 



i2 I f°° dA i\(s-s) 



1 ^ / _ n _ f/ .\2 / ..\2' 

'r'g dr'q exp < 

A<g<Ao ["A<|q|<Ao 

The rightmost integral, which is simply a product of separate Gaussian integrals, is equal to 



/ n '^r;<exp|l 0V-/3Kg^) + 



(18) 



n 



exp 



2 ^ In 

A<ij<Ao 



TT 



(19) 



^A<ij<Ao 

For large values of L, we may replace the discrete sum by an integral. Hence Eq. (19) can be 
approximated by 



exp 



J dq In yPnq — iXq' 



L r^o 

IT 



(20) 



up to a multiplicative constant independent of A, which we discard as it merely shifts F\R] by a 
fixed amount independent of R(x). We obtain 



exp{-/3F[R]} = exp -P ^ -k/|R, 



kl<A 



/oo - L f^o / \ 
dAexp iX{S — S) / dq In ( f3Kq^ — iXq"^] 
^oo vr 7a ^ ^ 



(21) 



yielding 



F[R] 



1 d^ny ^^^^ 

1 



dAexp iX{S — S) / dq\n\(3K.q^ — iXq 

-oo vr Ja ^ 



L r^o 



(22) 
(23) 



Thus the free-energy associated to R(x) contains a bending term [as previously for r(x)] plus a 
new potential (f>{S) function of the contour-length of the low-resolution filament's shape. Note that 
corrections to k would have been generated if we had taken into account the non-harmonic terms 
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in Eq. ( |T3| ) ( [Helfrich, 1985 ; Peliti and Leibler, 1985| ). A correct treatment of these terms, including 



measure factors ( Cai et al, 1994 ) is outside the scope of this paper. 

In the thermodynamic hmit (L oo), the integral in (p{S) can be approximated by its saddle- 



point value ( Walker, 1970 ). This yields 

iX{S) 



0(5) 



(S - 5) + A [""'dqln ((3Kq^ - iX{S)q^) , 



/5 



with X{S) the value of A that maximizes the integrand in Eq. (p3|); X{S) is defined by 

S-S _ 1 r^o dq 
L TT Ja p^q'^ - iX 



(24) 



(25) 



Note that X{S) actually lies on the imaginary axis of the complex plane. The Eqs. (|24| ) and (25) 
define (/>(5) implicitly. 
A quantity of interest is 

d4>{S) 



a{S) 



dS 



(26) 



It corresponds to the line tension associated with the low-resolution filament's shape R(x). Let us 
call (piS, A) the right-hand side of Eq. (^). Then we have 

~ dX 



d(t>{S,X{S)) ^d^ 
' dS dS dXdS' 



Since, by definition of A, d(j)/dX = 0, we obtain 



a{S) 



iXjS) 



(27) 



(28) 



Thus, the central result is the following: the line tension cr(5), the primitive of which is the potential 
(j){S), is given by 

S-S /-^o dq 



q^ H 

K 



(29) 



A. Analysis 

Let us introduce the following dimensionless quantities: 

S-S 



e 



tt/^kAq 
A 

a 



L 



kAI 



(30) 
(31) 
(32) 
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FIG. 2 Plot of the reduced effective wormlike chain's tension S against the reduced hidden length s. The 
hatched region is forbidden. The top curve corresponds to e = 10^^ and the lower one to e = 10^'^. For 
e < Smax there is no way to obtain a vanishing tension: the filament behaves as a tense "string" . 



The parameter s measures the hidden length of the filament. Since L < S < S, we have 

S-L 



< S < Sr, 



where Sr, 



ttPkAq- 



L 



(33) 



The dimensionless number PkAq, appearing in Eqs. ( pO| ) and (pq), is the ratio of the persistence 
length of the filament, Lp = (3k, to its thickness ~Aq ^. For a wormlike chain, we expect 



PkAo > 1. 



(34) 



The parameter e is the ratio of the high-resolution scale (the filament thickness) to the low- 
resolution scale A~^. Obviously, we shall consider the case 



e < 1. 



(35) 



The parameter S (having now a different meaning than in Sec. 1) measures the filament's tension 
in units of kAq. Note that for PkAq ^ 1, the quantity kAq is much larger than /cb^Ao, which we 
would expect if the tension were of molecular origin. 

In terms of these dimensionless quantities, the filament's tension, S(s), is defined by the equa- 
tion: 



dk 



fc2 + S' 



(36) 



10 



Let us now study the various regimes displayed by For s — > 0, E diverges as (1 — e)/s, as 

can be easily seen by neglecting k"^ with respect to S in Eq. (^). S decreases monotonically and 
vanishes at s = s*(e) with (for e <C 1) 

s* 0, (37) 

e 

it becomes then negative. When e < s^ax' ^i^^ beyond Smax and is thus unreachable (Fig.|^. 
Hence, as in the membrane case discussed by Fournier et al. (2001| ), there are two regimes: for 
e > s~axi i-6-> when the resolution A is higher than Ao/smaxj the filament is in a floppy regime, 
while for e < s~ax (low resolution), it is in a tense regime. 

In the floppy regime, the low-resolution contour R(x) exhibits a loose shape with an apparent 
contour-length fluctuating in the vicinity of the value S* (associated to s*) that minimizes the 
potential 0(5). At 5" = S*, we have a = 0: the tension effectively vanishes. This is reminescent of 
the case of membranes ( pavid and Leibler, 1991 ; Seifert, 1995| ). 



B. Effectively a tense string 

In the tense regime, the apparent contour-length minimizing 4>{S) corresponds to the lowest 
possible value of S, attained at s = Smax; which corresponds to the lowest possible value of S, i.e., 
L. The low-resolution, apparent shape of the filament is thus the straight /me joining the two fixed 
extremities (Fig.^. The filament, described at such low resolutions, behaves therefore effectively 
as a tense string of tension fio = c(L). As one pulls on the "string", the apparent contour- length 
increases {S > L) and so does the tension. 

Although Eq. ( p^ ) cannot be solved analytically for S, approximate solutions can be obtained 
in the tense regime both for medium and high tensions. Assuming e <C 1, we replace the lowest 
bound of the integral in Eq. (|36| ) by zero. This gives s ~ (l/\/S) arctan(l/\/S), which yields the 
following approximations: s ~ 7r/(2\/S) for < S <C 1 (medium tensions), and s ~ 1/S for S ^ 1 
(high tensions). 

For real wormlike chains and in the high tension case, we should actually take into account 
the finite stretching elasticity of the filament: this would modify our results. We shall therefore 
concentrate on the medium tension states of the tense regime. The relation s ~ 7r/(2-v/S) yields 

This represents the free-energy we have to pay in order to increase the apparent contour-length of 
the filament by a unit amount. It is worth noticing the intrinsic character of this expression: it 
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FIG. 3 The same wormlike chain viewed at high and low resolutions. The elasticity describing the low- 
resolution shape and its fluctuations (not drawn) is similar to that of a tense string. 

is independent both of the microscopic cutoff Aq and of the resolution A at which the filament is 
observed. For Eq. ( |38|) to hold, it is necessary that 

in order to be well into the tense regime (e ^ ■Sj^^x)- 
IV. CONCLUSION 

We have shown that at spatial resolutions A much less than the wavevector Ac defined above, a 
wormlike chain, i.e., a filament with a bending elasticity k and a fixed contour-length S, behaves 
as a tense string when it is attached at its extremities to a frame of length L. In other words, if one 
forgets about the microscopic wrinkles of the filament by considering only its low-resolution shape 
Yi[x) (the Fourier transform of the microscopic shape truncated at the cutoff A), one is actually 
considering a system that has the same type of elastic free-energy as a tense string. Indeed, in the 
minimum energy configuration, the low-resolution shape of the filament is straight and its tension 
is given by 

1 / L ^2 



which follows from Eq. (^) by setting S = L. Again, we note that this effective tension is inde- 
pendent of the resolution A (as long as A ^ Ac). For small deformations about this straight state, 
one can neglect the fact that a increases with S as one pulls on the "string". Then, the elastic 
free-energy of the filament, as seen at low-resolution, is 

F|R|.f.4.„(^)^ (41) 

In this last expression, we have neglected the remaining energy associated with the bending of 
the filament [see Eq. ([2^)]. It is actually negligible in the tense regime: for all the permitted 
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wavevectors {q < A), we have g <C Ac, and therefore aoq"^ » kq'^ since ^a^j k ~ Ac, as can be seen 
from Eqs. 



and (jig). 

Although this paper is focused on the concept of effective tension for a wormhke chain viewed 
at a low resolution, Eqs. (p2|), (p^) and (^) actually give the effective free-energy fit to describe 
an attached wormlike chain viewed at any resolution. 

It is interesting to compare the present results with those of previous works. For instance, using 
the wormlike chain model, Marko and Siggia (1995 ) have calculated the force / needed to extend 
a DNA molecule attached at its extremities. In agreement with experiments ( ^mith et al, 1992|) , 
they found that the end-to-end distance approaches the DNA contour-length as f~^^'^. Precisely, 
they obtained {Af Lp/k^T)^^/"^ ~ 1 — L/S in the strong force limit. Since /3k = Lp, Eq. ( pOD 
reproduces exactly the same result, upon identifying the effective tension o"o with the force /. 

Since the early work of Helfrich ( [Helfrich, 1985 ; Peliti and Leibler, 1985| ) , a large literature has 
been devoted to the renormalization of the elastic constants of membranes. No similar work has 
been performed on wormlike chains, however: this is in part because it is possible to directly calcu- 



late the tangent-tangent correlation function by analogy with the diffusion on a unit sphere ( Harris 
and Hearst, 196(^ ). A large amount of work has been devoted instead to the interpolation between 



the wormlike and the Gaussian chain ( Aragon and Pecora, 1985 ; Bawendi and Freed, 1985 ; Ha and 



Thirumalai, 1995 ). To the best of my knowledge, however, there has been no previous analysis in 
the spirit of the present paper. 
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